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Abstract 

We study in detail the dependence of the ^Sq superfluid pairing gap on nuclear 
interactions and on charge-independence breaking at the BCS level. Starting from 
chiral effective-field theory and conventional nucleon-nucleon (NN) interactions, we 
use the renormalization group to generate low- momentum interactions Mowfc with 
sharp and smooth regulators. The resulting BCS gaps are well constrained by the 
NN scattering phase shifts, and the cutoff dependence is very weak for sharp or 
sufficiently narrow smooth regulators with cutoffs A > 1.6 fm~^. It is therefore 
likely that the effect of three- nucleon interactions on ^Sq superfluidity is small at 
the BCS level. The charge dependence of nuclear interactions has a 10% effect on 
the pairing gap. 



1 Introduction 



Superfluidity plays a central role in strongly-interacting many-body systems 
ranging from nuclei, halo nuclei and neutron stars to cold atoms: The isospin 
dependence of nuclear pairing gaps shows striking trends over a range of iso- 
topes [1], the P decay of the two-neutron halo in ^^Li is suppressed due to 
pairing [2] similar to neutrino emission in neutron star cooling [3], and reso- 
nant Fermi gases exhibit vortices [1] and superfluid characteristics in thermo- 
dynamic [5] and spectroscopic properties [6]. 

For relative momenta k < 2fm~^, nucleon-nucleon (NN) interactions are well 
constrained by the existing scattering data [7]. The model dependence for 
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larger momenta shows up prominently, for instance, in the superfluid 
pairing gaps for Fermi momenta kp > 2fm~^ [8]. However, some uncertainty 
remains concerning a possible dependence of the ^Sq pairing gap on the in- 
put NN interaction in low-density neutron matter (kp < 1.6 fm~^). In this 
letter, we clarify this point and explore the dependence of ^Sq superfluidity 
on nuclear interactions at the BCS level in detail. We find that the BCS gap 
is well constrained by the NN phase shifts. Therefore, any uncertainties are 
due to polarization (induced interaction), dispersion and three-nucleon (3N) 
interaction effects. 

In addition to chiral effective-field theory (EFT) and conventional NN inter- 
actions, we use the renormalization group (RG) to evolve nuclear interactions 
to a lower resolution scale. The resulting class of low-momentum interactions 
^owA: [ZEUni, which is defined by a regulator with a variable cutoff A, re- 
produces the NN scattering phase shifts for momenta below A. We find that 
the cutoff dependence of the ^Sq BCS gap is very weak for sharp or suffi- 
ciently narrow smooth regulators with A > 1.6 fm~^. A comparison with the 
cutoff dependence found in A = 3,4 nuclei ("Tjon-line") [TT] and in nuclear 
matter [12j, when 3N interactions are neglected, suggests that 3N interaction 
effects on ^Sq superfluidity are small at the BCS level. 

Contact or separable pairing interactions can be implemented directly in 
current density-functional calculations. For low-momentum interactions, the 
weak-coupling approximation with a density-dependent contact interaction is 
reliable (see Ref. [13]), and a separable approximation is efficient [H]. There- 
fore, low-momentum interactions offer the possibility for a consistent treat- 
ment of the particle-hole and pairing channels in density-functional theory. 
Moreover, it is straightforward to adapt the RG to microscopically derive 
the renormalized pairing interaction introduced in the optimal regularization 
scheme of Bulgac [15] . 

This letter is organized as follows. In Sect. [H we present the formalism for 
^So superfluidity at the BCS level and our results for the dependence of the 
pairing gap on nuclear interactions and for the effects of charge-independence 
breaking. We conclude in Sect. [31 We emphasize that our work should be 
considered as a theoretical benchmark and not as a prediction of the superfluid 
pairing gap, since we do not include contributions beyond the BCS level. 



2 Formalism and results at the BCS level 

In the BCS approximation, the ^Sq superfluid gap A{k) is obtained by solving 
the gap equation with a free spectrum, e{p) = (in units c = h = m = 1, 



2 



with m the nucleon mass), 



where V{k, k') is the free NN interaction, ^(p) = e{p) — /i, and for a free spec- 
trum the chemical potential is given by /i = /c|/2. The Lippmann-Schwinger 
equation for the scattering T matrix in the same channel reads 

TT J L — 

Here and in the following, principal value integrals are implied. The homo- 
geneous gap equation can be understood as an equation for the residue of 
the pole of the T matrix, T{k,k';E) A{k)A{k')/[E - 2/x], for E 2fx, 
when the two-nucleon propagator in Eq. ([2]) is replaced by the corresponding 



self-consistent Nambu-Gorkov propagator 2/i — — * —2^^'^{j)) + /S?{p). This 
includes the propagation of back-to-back particle-particle, (1 — np)(l — rzp), 
and hole-hole modes, —UpUp, where Hp denotes the Fermi-Dirac distribution. 

For conventional large-cutoff and chiral EFT potentials, V{k, k') includes reg- 
ulating functions that render the integral convergent. These are of exponential 
form exp[-(A;2/A2)^] with A = 450 - 600 MeV in the current chiral EFT in- 
teractions at N^LO [T6p!7] . and phenomenological functions that imply large 
(few GeV) cutoffs in conventional NN potential models. 

For the RG evolution to low-momentum interactions Viowfc, we define a reduced 
interaction v{k,k') with V^owfc(^, = f{k)v{k,k')f{k'), where f{k) denotes 
a sharp or smooth regulator. Starting from an NN interaction with a large 
cutoff v{k, k') = VnnI^, k'), we use the RG equation [10], 
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(3) 



to generate low-momentum interactions with a variable cutoff A (both v and 
/ depend explicitly on A). Here, the reduced t matrix t(/c, /c'; E) is defined by 
Tiowk{k,k'; E) = f{k)t{k,k';E)f{k), where TjowA; is the solution to Eq. ([2]) 
with V = Viowfe- In this letter, we consider a sharp cutoff f{k) = 9{A — k) and 
smooth regulators of the exponential form f{k) = exp[— (/c^/A^)"], where n is 
a parameter that controls the smoothness. The resulting Viowfc is hermitian and 
preserves the low-momentum fully-on-shell T^n matrix, up to factors of the 
regulator function Tiowk{k, k; k"^) = p{k) T^^{k, k; k"^), as well as the deuteron 
binding energy [10] . 
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The RG equation, Eq. ([3]), is equivalent to a generalization of the Lee-Suzuki 
method p^flQ] and a subsequent Okubo hermitization [20] to smooth cut- 
offs [TUll^ . In the following, we will use the projection operator formalism to 
construct low-momentum interactions with a sharp cutoff and solve the RG 
equation for our results obtained with a smooth regulator. We note that the 
RG equation cannot be solved directly in the neutron- neutron ^Sq channel 
for most of the conventional NN interactions (except for the Nijmegen II po- 
tential), due to spurious resonances at high (~ GeV) momenta. The freedom 
in the choice of the regulator f{k) implies a scheme dependence of the gap 
A(/c) ~ f{k) at large momenta k ^ k-p (see Eq. ([T])). We will restrict our 
results to the gap on the Fermi surface A = A(A;f), where the momenta are 
on-shell and the gap is scheme independent. 

In the leading-order pionless EFT with sharp-cutoff regularization, one has 
V{k, k') = 9{A - k) elk - k')/[l/as - 2A/tt] (with scattering length a,). The 
resulting gap is cutoff independent for A ^ /i and large cutoffs, which follows 
from the gap equation with A{k) = 6{A — k) A, 



12. If, p 
A = ' J" 



A = --Jdp , (4) 



For A <ti fi and large A, the integral is given by 2 [—2kp + A + k-p ln(8/i/A)]. 
The UV divergence cancels against the cutoff dependence of the interaction 
in Eq. (jlj). This leads to the standard result A = 8fj,/e^ exp[7i/ {2kpas)] [22j . 

For the solution of the gap equation, we follow the method of Khodel et al. [23] : 
We first decompose the interaction into a separable and a non-separable part 

V{k, k') = V{kF, kp) (f){k) (j){k') + W{k, k') , (5) 

where (pik) = V{k,kF)/V{kF,kp) and W{k,k') is a suitably chosen function 
that vanishes when at least one argument is on the Fermi surface {k = kp). 
Then the gap equation, Eq. ([1]), can be replaced by an equivalent system of 
two equations. 



m=xik)+'-[dpp^^M^, (6) 



/eiP) + A'x'ip) 

= 1 + V{kp,kp) - dpp (7) 



e{p) + A^x'{p) 

where A{k) = Ax{k) with x(^f) = 1- This system has the advantage that 
the integrand in Eq. ([6]) vanishes on the Fermi surface, and consequently the 
function x(^) is only weakly sensitive to changes of A{p) in the denomina- 
tor. Therefore, to a good approximation, Eq. can be linearized. In the 
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Fig. 1. The neutron- neutron Sq superfluid pairing gap on the Fermi surface 
A = A(A;f) versus Fermi momentum kp for low-momentum interactions Vjowfc with 
a sharp cutoff A = 2.1 fm~^. Mowfc is derived from various charge-dependent NN 
interactions |16|24|25|26j . We have verified that the results are cutoff independent 
from A = 1.6 fm~^ to A = 2.5 fm~^. The inset magnifies the smah dependence on 
nuclear interactions near the maximum. 

first iteration, we solve Eq. ([6]) by inversion using a sufficiently small constant 
Ansatz for A{p) in the denominator. Next, we solve Eq. ([7]) for A using New- 
ton's method with the solution for x(^) of previous step. The iteration 
of this procedure (where A(p) in the denominator of Eq. is updated at 
each step) leads to a rapidly converging solution for the BCS gaps. We follow 
this method for the smooth-cutoff V^owfc- For the sharp-cutoff case, we have 
found that a solution of the gap equation, Eq. ([1]), by iteration with an initial 



guess A{k) ~ \j\kFV\owk{k,k)\ also converges rapidly. This Ansatz is moti- 
vated by the observation that a separable approximation is meaningful for 
low-momentum interactions [14j. 

Our results for the density dependence of the neutron-neutron ^Sq superfluid 
gap A are shown in Fig. [H The low- momentum interactions Viowfc are derived 
from various charge-dependent NN potentials [T6ll24f25ll26] using a sharp cutoff 
A = 2.1fm-\ We find that the BCS gap is almost independent of the NN 
interaction. Consequently, we conclude that the ^Sq gap is strongly constrained 
by the NN scattering phase shifts. This has been noted previously (see for 
example Ref. [27j), but without considering charge dependences. Moreover, 
these are the first results for chiral interactions at N^LO. We use the N'^LO 
chiral potential of Ref. fL6\ , since it is the chiral interaction that leads to the 
most accuracte reproduction of the phase shifts. 

The maximal gap at the BCS level is A ^ 2.9-3.0 MeV for kp ^ 0.8-0.9 fm"^ 
The small deviation of the N^LO gap from the band at higher densities in 
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Fig. 2. The charge dependence of the Sq superfluid pairing gap A versus kp. The 
Unes indicated in the legend are the neutron-proton gaps, whereas the grey lines 
show the neutron-neutron gaps from Fig. [TJ For further details, see the caption of 
Fig. [TJ We have also verified that the neutron-proton gaps are cutoff independent 
over the same sharp-cutoff range. 

Fig. [1] is consistent with the slightly more attractive ^Sq phase shifts at the 
corresponding energies (compare, for example, the phase shifts of the CD- 
Bonn [26] and N^LO potentials). We find that the gaps are cutoff independent 
over the range considered here, A = 1.6 fm~^ to A = 2.5 fm~^. This result 
is consistent with the findings of Kaiser et al. that the cutoff dependence is 
substantially reduced for chiral EFT interactions when going from NLO to 
N^LO [28j, since the latter leads to a better description of the NN scattering 
phase shifts. In addition, the BCS gaps for the "bare" interactions are within 
2% of the V\awk results shown in Fig. [T]for k-p < 1.0 fm~^, and the difference 
is compatible with the spread in the V^iowfc result over all densities. (This also 
holds for Fig. [21) For completeness, we mention that Sedrakian et al. [29] have 
solved the BCS gap equation for one low-momentum interaction (Viowfc derived 
from Nijmegen 93 ^24j with A = 2.5 fm~ ), but they did not explore the cutoff 
dependence. 

Isospin symmetry breaking leads to small charge dependences in nuclear inter- 
actions. As a result, the ^Sq neutron-proton scattering length Onp = —23. 768 ± 
0.006 fm [3^ is more attractive than the neutron-neutron scattering length 
= —18.5 ± 0.3 fm [21] in the same channel. This effect is dominantly due 
to the charge dependence of the one pion-exchange interaction The cen- 
tral part of VJr in the neutron-proton charge-exchange channel is of the form 
— m^±/(g'^ -|- m^±), where q' is the (exchange) momentum transfer. Since 
the charged pion is heavier than the neutral one, m^± = 139.57 MeV and 
m^o = 134.98 MeV, the resulting neutron-proton interaction is more attrac- 
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Fig. 3. The neutron-neutron So superfluid pairing gap A as a function of the cutoff 
A for three densities and different smooth exponential regulators, as well as for a 
sharp cutoff. The low-momentum interactions are derived from the N^LO chiral 
potential of Ref . [TB] . 
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Fig. 4. The neutron-neutron Sq superfluid pairing gap A as a function of the 
cutoff A for three densities. The low- momentum interactions are derived from the 
N^LO and the Nijmegen II [24J potential with exponential regulator n = 7. 

five. In Fig. [21 we show the charge dependence of the ^Sq superfluid pairing 
gap versus Fermi momentum. The neutron-proton gaps are ~ 0.3 MeV larger 
at maximum with a slight shift to higher densities. The 10% effect on the 
pairing gaps clearly reflects the charge dependence of nuclear interactions. 

Next, we study the dependence of the neutron-neutron ^Sq superfluid pairing 
gap as a function of the cutoff starting from the N^LO chiral interaction. 
Our results for three representative densities and different smooth exponential 
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regulators f{k) = exp[— (fc^/A^)"], as well as for a sharp cutoff, are shown in 
Fig. [3l As long as the cutoff is large compared to the dominant momentum 
components in the Cooper bound state, the gap depends very weakly on the 
cutoff. Below this scale, which depends on the density and the smoothness 
of the regulator, the strength of the bound state decreases, since some of 
the momentum modes that build up the Cooper pairs are integrated out. 
From Fig. [3], we observe that the cutoff dependence is very weak for sharp or 
sufficiently narrow smooth regulators with A > 1.6 fm~^. It can be seen that 
n = 3 is too smooth, but that n > 5 is sufficient. For lower densities, even 
lower cutoffs with A > 1.2kp are possible. 

The N^LO chiral interaction has a cutoff A = 2.5 fm^^ (or 500 MeV) ^6] and 
one may suspect that the cutoff dependence could be larger for conventional 
NN potentials. In Fig. |H we show that this is not the case by comparing 
the gaps from Fig. [3] to results obtained with the Nijmegen II potential |24j . 
which has a large (~ GeV) cutoff. The resulting cutoff dependences are similar 
and in particular very weak for sharp or sufficiently narrow smooth regulators 
with A > 1.6 fm~^. This shows that the ^Sq superfiuid pairing gap probes 
low-momentum physics. 



3 Conclusions 

In this Letter, we have presented a systematic study of the interaction depen- 
dence of the BCS ^So superfiuid pairing gap. We have shown that this gap 
is practically independent of the choice of the NN interaction, and therefore 
well constrained by the NN phase shifts. Furthermore, we have found only 
a very weak dependence on the cutoff for low- momentum interactions Viowfc 
with sharp or sufficiently narrow smooth regulators for A > 1.6 fm~ . For low 
densities, it is possible to lower the cutoff further to A > 1.2A;f- We also find 
that the pairing gap clearly reflects the charge dependence of NN interactions. 
Neutron- neutron and neutron-proton ^Sq are not carefully distinguished in pre- 
vious work. We conclude that the uncertainties in ^Sq superfluidity are due 
to an approximate treatment of induced interactions and dispersion effects, 
which go beyond the BCS level, as well as due to 3N interactions. 

The weak cutoff dependence indicates that, in the ^Sq channel, the contribu- 
tion of 3N interactions is small at the BCS level. We note that the 27r-exchange 
part ( "Cj-terms" ) of the corresponding low-momentum 3N interactions [TT] is 
cutoff independent up to the regulator functions. The latter lead to a cut- 
off dependence as the density increases (see, for example, the Hartree-Fock 
results in Ref. [12]). We also emphasize that 3N interactions will contribute 
differently to ^Sq superfluidity in pure neutron matter compared to symmet- 
ric nuclear matter. Additional insights will come from an investigation of the 
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cutoff dependence of ^Si-^Di superfluidity, wliere 3N interactions may play 
an important role. Work in this direction is in progress [32] . 

Low-momentum interactions, via weak-coupling or separable approximations, 
can be implemented directly in current density-functional calculations. Fur- 
thermore, it is straightforward to adapt the RG used here to a microscopic 
derivation of the optimal pairing interaction of Ref. |15j . 

Finally, we emphasize that our results are obtained at the BCS level and do 
not include polarization and self-energy effects. Therefore, our work should be 
considered as a theoretical benchmark and not as a prediction of the super- 
fluid pairing gap. Recently, polarization effects on the pairing gap have been 
studied in an RG approach to the many-body problem, starting from low- 
momentum NN interactions [13]. Is is found that polarization effects lead to 
a suppression to a maximal gap A f» 0.8 MeV, in qualitative agreement with 
the earlier work of Wambach et al. [33] . The RG approach is nonperturbative 
and includes long-range particle-hole induced interactions and the dominant 
self-energy effects. Tensor induced interactions are not included in Ref. [13] 
because their effects on ^Sq superfluidity are expected to be small at very low 
densities, and because the strength of the tensor force is coupled to the cor- 
responding 3N interaction [TU|TT] . In addition, there are recent finite-particle- 
number, fixed-node AFD Monte Carlo calculations [51], which in principle 
include olarization effects, but result in a maximal gap A k, 2.5 MeV. How- 
ever, in these calculations larger particle numbers may be required to capture 
long-range polarization effects. 
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